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r"! ■ Abstract: A simple recursive expansion algorithm for the integrals of superstring five point 

amplitudes in a flat background is given which reduces the expansion to simple symbol(ic) 
CN ! manipulations. This approach can be used for instance to prove the expansion is maximally 

transcendental to all orders and to verify several conjectures made in recent literature to high 
order. For closed strings in particular the maximal R-symmetry violating amplitudes in type 
IIB superstring theory are studied. The obtained expansion of the open string amplitudes 
reduces the analysis for MRV amplitudes to the classification of completely symmetric poly- 



C^ ■ nomials of the external legs, up to momentum conservation. Using Molien's theorem as a 



counting tool this problem is solved by constructing an explicit nine element basis for this 
class. This theorem may be of wider interest: as is illustrated at higher points it can be used 
to calculate dimensions of polynomials of external momenta invariant under any finite group 
for in principle any number of legs, up to momentum conservation. We comment briefly on 
extension of the techniques discussed in this paper to higher numbers of legs. 
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1 Introduction 

The ability to calculate is central to achieving understanding in theoretical high energy 
physics. For the specific case of scattering amplitudes, progress can be measured in gen- 
eral along the twin axes of numbers of loops and legs. In string theory, a third axis appears: 
that of orders in the field theory expansion. Whenever the momenta in a particular scattering 
experiment are small compared to the string scale, a perturbation theory may be set up. In 
practice, this is often referred to as the ol expansion. Furthermore, it is part of string theory 
textbooks to calculate scattering amplitudes in a flat background at the string tree level in 
the form of integrals over the boundary of the string worldsheet^. In principle the field theory 
expansion of the amplitudes follows from the integrals. Beyond the four point case however 
relatively little is known about these string scattering amplitudes and especially their field 
theory expansion. This is quite surprising, since these scattering amplitudes are one of the 
few things we know how to compute reliably in string theory and contains much information 
on the theory. This paper is part of a wider effort to address this unsatisfactory situation. 

Considering the birthplace of string theory [4] it is humbling to think that the first 
systematic attempt at computing a five point superstring amplitude appeared only relatively 



^In the closed (or mixed) case this follows after apphcation of the Kawai-LeweUen-Tye [1] relations (or their 
analogons [2], [3]) 



recently in [5]. A first all multiplicity result for the first string theory correction at order 
(a')^ appeared in [6] in four dimensional kinematics. More simplifications for the integrand 
of the amplitudes were obtained in [7], where a particularly efficient way of ordering the 
open string theory amplitudes in terms of their field theory limits was introduced. Some first 
results for the structure of the field theory expansion appeared in e.g. [8], [9] and [10]. A first 
systematic study of the field theory expansion beyond simple orders however was announced 
only very recently [11]. The structure found there 'experimentally' is currently only partially 
understood, see e.g. [12]. One goal of this article is to reproduce these results and push the 
investigation to higher orders in ol . 

It should be noted that it is known that for five point amplitude the functions which ap- 
pear in the superstring are of the 3-F2 hypergeometric type. Expansion of this type of function 
around integer arguments is a classical problem which appears often during the evaluation 
of complicated Feynman integrals in the more phenomenologically oriented sections of high 
energy physics. The literature on this is correspondingly vast: see for instance the references 
in [13] and [14]. The papers [13] and [14] in particular introduce an expansion algorithm for 
hypergeometric functions which will be implemented in this article for superstring amplitudes 
with five massless legs. This algorithm highlights the close connection of the expansion to 
certain poly logarithms. This is seen as a good starting point to deepen the investigation of 
the field theory expansion of the superstring. The techniques covered here do, in principle, 
generalize to higher leg multiplicity. Note that a naive field theory expansion does not make 
sense for amplitudes which involve any of the Regge excitations: their masses are given in 
string units l/(a'). 

This paper is structured as follows: In section 2 some notation is established by writing 
several recently made conjectures. The expansion algorithm for open string amplitudes is 
subject of section 3. In section 4 the focus is on the closed string amplitudes which follow 
from the KLT relations. Special attention is focused on a recently identified class of partic- 
ularly simple amplitudes: the maximal R-symmetry violating ones. It is shown these take 
in the case at hand a particularly simple expression in terms of the open string amplitudes. 
To facilitate the discussion a simple nine element basis is introduced for the completely sym- 
metric polynomials of external momenta. The used technology quite likely has applications 
beyond the immediate application at hand, as will be highlighted at the end of the section. A 
discussion and conclusion section round of the main presentation. The main results of this pa- 
per (explicit matrices, expansion coefficients, (/)-map) are contained in the archive submission 
for the readers convenience. 

Note aided in "proof 
While this paper was being readied for publication an unexpected overlap was discovered with 
a project by different authors which also aims at novel expansion methods for superstring 
amplitudes. Results from this project have since been announced in [15] and [16]. 



2 Brief review of five point superstring ampfitudes 

Many of the conjectures in this section originate in [11]. Superstring amplitudes are functions 
of the external kinematics and polarizations. As shown in [7] the polarization dependence of 
open superstring amplitudes resides completely in Yang-Mills tree amplitudes. The rest is a 
relatively simple function of the momenta of all particles k^ subject to momentum conserva- 
tion 

n 

E^f = (2-1) 






for n incoming on-shell and massless (a'(fcj)^ = -gSa = 0) momenta. Lorentz invariance then 
dictates that the basic variables are made of 

Sij = a' {h + kjf (2.2) 

where a' is the usual dimension-full parameter which set the string scale. The dimensionless 
Sjj-variables will be referred to as Mandelstam variables or 'Mandelstams' for any number 
of external legs. Momentum conservation implies relations between the Sij. Let us first 
assume that the number of dimensions is larger than or equal to n. Then constraints on the 
Mandelstams can be derived from momentum conservation as displayed in equation (2.1) by 
contracting with all n possible momenta. Hence there are 

#vars = -n{n — 3) (2-3) 

independent variables. For four points, these boil down to Mandelstam's original variables, 
subject to s + t + n = 0. For five points an interesting set for our purposes is the cyclicly 
symmetric set 

{Sl2, S23, 534, S45, S51} = {si, S2, S3, 54, S5} (2.4) 

The other 5 Mandelstams can be expressed in terms of these by solving the momentum 
conserving equations as 



(2.5) 



For general multiplicity n a set of independent variables can be chosen to be for instance 

{sij \i < i < n - 2} U {si,„_i I 1 < i < n - 3} (2.6) 

which corresponds to the dependent variables 

{Si,n\i<n-l} U {Sn-2,n-l} (2.7) 



'Sl3 = 


-S12 


- S23 + 545 


Sl4 = 


-S5I 


+ S23 


- S45 


S24 = 


S5I - 


S23 - 


S34 


S25 = 


-S12 


- S5I 


+ S34 


. S35 = 


S12 - 


S34 - 


S45 



It is straightforward to check that the latter set of variables has a solution in terms of the 
former set under momentum conservation. Other sets are possible of course, but will not be 
needed here. There are more constraints when the number of legs is more than the number of 
dimensions plus one. This is simply due to the fact that n vectors in a D dimensional space 
are always dependent when n > D. Momentum conservation accounts for the '+1' offset. 
Since the subject of this note are superstring amplitudes, this will discrepancy will only kick 
in for 11 or more particles and can safely be ignored for the present purposes. 

2.1 Open superstring 

The two independent five point open superstring amplitudes from which the others can be 
derived can be chosen to be the color-ordered amplitudes A°(12345) and j4°(13245). These 
can be expressed in terms of the Yang- Mills amplitudes with the same ordering as 



.4°(12345)\ ^ /ylYM(i2345)\ 

A° {13245) J yA^^ (13245) J ^ ' 



where 




(2.9) 



with the entries defined as the integrals 

-1 rl 

'0 Jo 



Fi = S12S34 / / dxdy x'^'^y'^^'^il - x)"34-i(i _ y)^23(i _ ^y)^24 (2.IO) 

Jo Jo 



and 



-^2 = 513524/ [ dxdy x'^"y''%l-x)''*{l-y)'^^{l-xy)'^*-^ (2.11) 

Jo Jo 

and the functions Fi are the Fi with particles 2 and 3 interchanged. 

MZVs and polylogarithms 

The expansion of these functions in F or indeed their higher point generalizations contains 
generically many Multi-Zeta Values (MZVs). These are certain iterated integrals, specified 
by a vector of "depth" d whose entries are positive integers that sum to the so-called "weight" 
w. 

C(ni,n2,...nrf) (2.12) 

d 
^nj=u; (2.13) 

i=i 

Given a vector n one can construct a ordered vector n with entries or 1 such that there are 
rii — 1 consecutive O's follows by a one. 

n = {0,0,...,l,0,...,l,...} (2.14) 



The multi-zeta values are very well studied and obey many relations with values in the 
rationals Q. An explicit basis up to weight 22 exists for these [17]: this basis will be referred 
to as the datamine basis. Of central importance in this article is the realization that the 
multi-zeta values are boundary values of polylogarithms. To write these define an iterated 
integral as 

I{a):f{x)^rdt-^f{t) (2.15) 

Jo *-a 

which takes a function f{x) and integrates it to a new function g{z). In this notation, the 

multi-zeta value is defined by repeated application of the integral operators, specified by the 

vector n, 

C(fi) = [/(Or-il(l)/(0)"^"i/(l) . . . 1]^^, (2.16) 

while the boundary value of the last integral is set to 1. In more modern language, the vector 
n is the symbol of a particular harmonic polylogarithm. The multi-zeta values appear as 
'boundary values' of the polylogarithms. In this context it should be noted that the integrals 
appearing in string amplitudes can always be expressed as boundary values of functions as 
noted long ago [18] . 

Definitions, conjectures 

Any known expansion algorithm for the integrals appearing in the five point amplitude will 
yield the answer in terms of generic MZVs. Using the datamine basis^ then gives an unam- 
biguous answer which may be compared between approaches. 

In terms of this basis one can now define two by two matrices Pi , Mj as the coefficient of 

P2i = -^lc(2j){a')2^ (2-17) 

and 

M2i-1 = -^|^{2i-l)(a')^'-^ (2.18) 

Conjecture: in terms of the a' expansion the function F factorizes as 



F = B[Pi]C[Mi] conjecture (2.19) 



with 

B = (f2i^TC2P2i] (2.20) 



.j=0 



Moreover, coefficients exist such that 



C{Mi) = ^(a)'' ^/j(M(^''^')) conjecture (2.21) 



^Not forgetting that two MZVs with fii and h2 are the same whenever one yields the other when inter- 
changing ones and zeros and inverting the order. This is known as the 'duahty' relation. 



where M^^'^' ranges over all products of Mi matrices such that in the product 

\{Mi,-...Mi. (2.22) 

3 

the restriction that 

^ij = k (2.23) 

i 
holds, with the convention that Mcven and Mi are zero. A second conjecture states that the 
coefficients / follows from a map from the multi-zeta values to non-commutative polynomials 
given by the map constructed in [19]. The non-commutative polynomials are constructed 
out of non-commutative products of monomials /2i+i for i > and the commutative element 
/2 ^ C,2- The conjecture is that for a certain product of M matrices the coefficient is the 
noncommutative polynomial with the same indices as the M's but the order inverted. The 
required map map can at least in principle be constructed recursively (see [19] for details), 
although in practice a basis such as the datamine one is very useful to cut down on the needed 
algebra. Results for the (p map up to weight 21 are included in the archive submission of this 
paper. Taken together these conjectures, if true, imply that once the matrices M and P are 
known the full superstring amplitude can be reconstructed. 

2.2 Closed superstring 

The KLT relations generically express closed superstring amplitudes in terms of open super- 
string amplitudes at the tree level. In the five point case this can be conveniently written 
as 

,_/^°(12345)\^ /a°(12345)\ 



where 



Cll C\2 
C2I ^22 



(2.25) 



with coefficients given in [11]. The field theory limit of this matrix will be denoted Sq. Two 
conjectures were made in [11] about these: 





P^SP = So 


conjecture 


(2.26) 


and 




Mf So = SoAfi 


conjecture 


(2.27) 



The conjecture in equation (2.26) was verified up to and including order 17, while the conjec- 
ture in equation (2.27) was verified to and including order 16. The first equation especially 
yields a significant simplification of closed string amplitudes. 



Fi = 



3 Field theory expansion of open string five point ampHtudes 

In this section we show how the conjectures recalled above can be checked to high order. 
The driver here is the ability to expand the 3F2 hypergeonietric function which appear in 
the evaluation of the integrals in equations (2.10) and (2.10). A neat way to obtain the 
hypergeonietric function out of the integral is to use a Mellin-Barnes representation for the 
"(1 — xy)^" type term in the integrand which reduces the integral tobeta function type 
integrals. Performing this integral and re-writing the result as an infinite sum shows the 
coefficients of the integrals are such that the sum is a hypergeometric 3^2 function. In 
particular, 

, ^ r(i + si)r(i + s2)r(i + s3)r(i + s4) j. / ^i, 54 + 1,^2 + 53-^5 . \ , .3^^ 

^ r(l + Si + S2)r(l + S3 + S4) ^ ^ \1 + S1+S2,1 + S3 + S4' 7 ^=^ 

p f w , r(i + gi)r(i + g2)r(i + g3)r(i + s4) 

F2 = [Si -Si- S2){S5 -S2- S3) — — ^ 

r(2 + si + S2)r(2 + S3 + S4) 

M +Si,l + S4,l + S2 + S3 -S5 V 

\ 2 + si + S2,2 + S3 + S4 / 

hold. Note that the expansion of the pre-factor in both cases is closely related to the four 
point amplitudes. These functions turn out to be a nice warmup example for expanding the 
3F2 functions. 

The technique for expanding hypergeometric functions which appears in this section was 
first proposed in the field theory context in [13] and [14]. The reader is referred to these papers 
for a full explanation. In the course of this project also the HypExp 2.0 [20], XSummer [21] 
and NestedSums [22] packages written in Mathematica, Form and C++ respectively were 
evaluated for the purpose of expanding the hypergeometric functions in equations (3.1) and 
(3.2). For this particular application however these were slower by many orders of magnitude 
than the algorithm explained below. 

3.1 Warmup: expanding the 2F1 hypergeonietric function 

Of interest here is the four point amplitude which follows from the 2-^1 hypergeometric func- 
tion as: 

/, / , / M r(aa' + l)r(6a' + 1) ,„ „, 

2F1 -aa', 6a'; 1 + ba'; z) ,=1 = \, ,'\^,. ' 3.3 

Y[aa' + ba' + 1) 

where a' dependence has been emphasized. Hypergeometric functions generically obey dif- 
ferential equations. In this case this equation reads: 

{z{e - a' a){e + a'b)- 9(9 + a' b)) 2F1 (-aa', ba'; 1 + ba'; z) = (3.4) 



where 9 = z-^. As above, we are eventually interested in the boundary value z = 1. The 
trick is now to first expand in a' before taking a' -^ 0. Write for this expansion 

oo 

2Fi{-aa'M';l + ha';z) = ^Wk{z){a'f (3.5) 

fc=0 

With this expansion the above differential equation can be written as a recursive differential 
equation, 

{z - l)e'^Wk - {z{a -b) + b) Owk-i - zahwk~2 = (3.6) 

The first few terms can be determined directly from the expansion (or the integral): wq = 
l,w;i = 0. Furthermore, there are boundary conditions: 

Wk{<d) = Q {9wk){<d) = Q yk>0 (3.7) 

With these boundary conditions the above equation can be solved by integrating the equation 
using the integrands suggested by the harmonic polylogarithms. This reduces to a 'symbol- 
level' recursion relation for the expansion of 2-^1 (—aa', ba'; 1 + ba'; z): 

Wj = -a({0, 1} U {wj^i)\l) + 6({0} U Wj^i) - a6({0, 1} U Wj-2) (3.8) 

Where the join operator is acting linear on the symbols of the polylogs appearing in the 
previous coefficient and the \1 is a linear map which drop the leading entry of the symbol. 
This is the action of 9 on an integral with symbol entry 0. From this recursion relation one 
obtains immediately that the expansion of this function is maximally transcendental: the 
multi-zeta values which appear at order k in the expansion have weight k if those at weight 
k — 1 do. Calculating the field theory expansion of the four point amplitude at a particular 
order now boils down to implementing the recursion relation and setting z = 1 after the 
last recursive step. Even in Mathematica, this can be made to yield high orders in very 
little time (order 100 needs about 1.5 sec on an ordinary laptop). Furthermore, some simple 
experimentation yields the result that the coefficient of a given a^b^ term in the expansion 
has the intriguing form 

~aVC({0*,F}) (3.9) 

with unit coefficient. On other grounds, see e.g. [11], it is known that the expansion of the 
four point amplitude in the string theory can be expressed in terms of single zeta values only. 
This suggest a theorem that the above MZVs can be expressed in terms of single zetas only. 
It turns out that this theorem has appeared in the math literature in [23]. 

3.2 Expanding 3F2 hypergeometric functions 

Now the same strategy can be employed for the 3F2 hypergeometric functions in equation 
(3.1) as well as in equation (3.2). Fundamental here is the expansion of the function 

3F2[ ,, ,, ■z\='^wk{z)a (3.10) 



a'bi, a'b2 



The needed hypergeometric functions in the five-point amphtudes are related to this one by 
the relations 

I 1 + a'ai,l + a'a2,l + a'as ^ \ _ (1 + a'6i)(l + a'62)^ p I a'ai,a'a2,a'a3 _ \ 
I 2 + a'6i, 2 + a'62 ' I {a')'^aia2a3Z I a'h\,a'h2 ' / 

(3.11) 
and 

I l + a'ai,a'a2,a'a3 \ 1 //i , / \ ^ ( a'ai,a'a2,a'a3 \ 

\ 1 + a oi, 1 + a 62 / a oi \ a 61, a 62 / 

In passing we note that these relations show that the two by two matrix F can be schematically 
be factorized as 

/ r(i+s2)r(i+s3)r(i+s4)r(i+s5) ^ \ 
p ^ r(i+s2+s3)r(i+s4+s5) 

V (2^3)^ 

((3F2 \ /^7^^^3i^2^7^^^3i^2\\ 

1^1^ (2o3)y ^ (2o3) (203)^^''' 

where the bottom line is the result of the exchange of particles 2 and 3 of the top line, with 
the line reversed. 

To write the analog of the recursion relation of equation (3.8) above it is (also programming- 
wise) useful to introduce new variables for the function in equation (3.10) 

Ai = ai + 02 + 03 - 61 - 62 A2 = aia2 + a2a3 + «3ai - 6162 A3 = 010203 

Qi = bi + &2 Q2 = hh 

which make obvious the usual symmetries of the function. Inverting the differential equations 
with boundary conditions then yields the recursion relation 

Wj = Ai({0,0,l} U {wj.i)\2) + A2({0,0, 1} U (u;,-2)\l) 

+ A3({0, 0, 1} U (wj^s) - Qim U (^,„i) - Q2({0, 0} U (u;,_2) (3.14) 

with wq = 1 and wi = W2 = a,s starting conditions. Recalling that the action of the 9 
operator simply removes a leading from the symbol of the polylog then yields a formula 
to expand the hypergeometric functions in equations (3.1) and (3.2). Note that this makes 
it clear that the resulting expansion is maximally transcendental: if the expansion has this 
property to order k, then the above relation extends this to k + 1. The base step in this 
inductive argument is trivial. A second feature of the above expansion is that since the fourth 
coefficient is proportional to A3, every coefficient in the expansion is. This can be used to 
verify that the expansion of the function in (2.10) and (2.11) is really a strict polynomial. 
Combining the recursion relation for the 3-F2 hypergeometric function and its pre-factor is an 
interesting problem left to further research. 



3.3 Verifying conjectures up to and including weight 21. 

Calculating the field theory expansion of the five point amplitude at a particular order now 
boils down to implementing the recursion relation and setting z = 1 after the last recursive 
step and after taking any required z derivatives. There is a FORM implementation of this 
algorithm [24] which has been used to get the expansion up to and including order 28, in 
terms of series of terms of quite general MZVs. A less quick but still sufficient Mathematica 
implementation was used instead to verify the conjectures made by [11], recalled above, up 
to and including order 21. It should be said that the bottleneck here is not expanding the 
hypergeometric functions but instead is the substitution of the datamine basis for the MZVs 
which creates memory issues'^. 

In general the landscape of results does not change from that sketched up to and including 
order 16 [11]. New is the appearance of relations amongst products of M matrices, first at 
order 18. These look remarkably like Jacobi-like relations for some algebra, but we have been 
unable to identify this algebra. At order 18 the one relation appearing takes the form: 

- M3 • M5 • M3 • M7 + M3 • M5 • M7 ■ MS + M3 ■ M7 ■ M3 ■ Mb 

-M3-M7- Mb ■ M3 + Mb ■ M3 ■ M3 ■ M7 - Mb ■ M3 ■ M7 ■ M3 

- M7 ■ M3 ■ M3 ■ Mb + M7 ■ M3 ■ Mb ■ M3 = (3.15) 

This relation was elegantly explained in [12]. Similarly, there is one such relation at order 19, 
three at order 20 and six at order 21. This lead to a further practical remark: in checking 
the conjecture of equation (2.21) it turns out to be handy to first take the trace and then 
write the conjectured sum on the right hand side as a vector with the trace of each matrix a 
separate entry. Then this vector can be evaluated on the same number of data points as the 
length of the vector. This reduces inverting the problem to inverting a matrix, which is easily 
implemented in Mathematical. Care has to be taken though to verify that not an accidentally 
degenerate set of data points was selected. Explicit results for the expansion coefficients as 
well as the M and P matrices are included in the archive version of this submission. 

One discrepancy remains at order 21 after verifying the conjectures recalled in equations 
(2. 19), (2. 21), (2. 26) and (2.27). This has to do with Brown's (jy-map which is supposed to 
map all coefficients from the datamine-basis to simple monomials in terms of the /2i+i type 
variables. Even taking into account the relations between products of M-type matrices at 
order 21 we have been unable to obtain this result for those matrix combinations which contain 
one Mg and 4 M3 matrices and a class with M^jMj and M3. Details can be reconstructed 
out of the supplied data in the archive version of this submission, or by asking the author. 
This might have it's root in the exact procedure used to regulate some of the integrals (with 



^This restricts to roughly order 16 on a machine with 8GB of memory for a calculation which takes about 
10 minutes in Mathematica, much of which is spanned compiling MZV data. 

* "LinearSolve" is your friend here. Note this friend can solve equations involving non-invertible matrices 
(if the system has a solution). For speed, invert for every MZV appearing separately. 
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leading 'ones') which appear in Brown's procedure^ or less high-brow in a calculational error 
not yet identified. It will be interesting to locate the source of this discrepancy. 

4 Field theory expansion of closed string five point amplitudes 

By the KLT relations, the previous results on field theory expanding the open superstring 
amplitudes have an immediate application in the field theory expansion of closed string am- 
plitudes. Especially conjecture (2.26) makes clear that there are large cancellations hidden 
in the KLT relations. Although true, the resulting expressions are generically complicated 
and not very illuminating without further insight. In this article the closed string amplitudes 
are therefore considered in a special sub sector within the type IIB string in ten dimensions 
uncovered in [25] which is remarkably similar to MHV amplitudes in four dimensions. This 
class of amplitudes violates R-symmetry maximally, whereas the four dimensional MHV am- 
plitudes violate helicity maximally. Structurally, these amplitudes can be expressed on an 
appropriate [26] on-shell superspace as 

where j4 is a completely symmetric function of the external momenta which does not have 
massless poles for more than four particles. This implies in the field theory expansion that 
this function has no poles: it must be a strict polynomial. 

By Bose symmetry MRV amplitudes in the field theory expansion are therefore expressible 
in terms of completely symmetric polynomia of the external momenta, as the fermionic delta 
function is completely symmetric. Lorentz invariance dictates that these must be polynomials 
in the Mandelstams, Sjj. The problem now is momentum conservation: solving this explicitly 
yields variables which generally do no transform nicely under permutations. This is already 
true in the four point case: one element of the set s,t is, under permutations of the external 
legs certainly mapped to u. For four particles it is known that any completely symmetric 
polynomial / can be written as a double polynomial expansion in two basis elements, 

f{s,t,u) = ^c,jaiai (4.2) 

where 

(J2 = s^ + f + u^ 0-3 = (stu) (4.3) 

This will be re-derived below. 



^Technically, in the current implementation some symbols appearing in the map with 'leading ones' which 
are divergent are reduced to leading zero type by using, in Brown's notation, /(O; 1; 1)=0, and the fact that 
the integrals satisfy a shuffle algebra. This regularization was suggested in [11] 



11 



4.1 Classifying completely symmetric polynomials 

To the best of our knowledge a similar basis is unknown above four points in the existing 
literature. One can, in principle, construct completely symmetric polynomials by brute force 
by summing certain base polynomials over all permutations. Relations between these could be 
established by evaluating these, say, k polynomials at k different kinematic points by taking 
Sij to be natural numbers. The resulting k hy k matrix has null eigenvectors which are the 
sought-for relations. Implementing this strategy for five points and comparing to the online 
integer sequence database [27] then led to a far better approach through Molien's theorem. 
Molien's theorem might be of wider interest so is presented here in a general formulation: 



Theorem 4.1 (Molien) Let V be a vector space and M a matrix representation of a group 
P, which is a subgroup of the permutation group on this space. Let fk be a function on the 
space V invariant under the action of M , with fixed homogeneity e.g. 

fk{\v) = X'fkiv) (4.4) 

Let dk be the number of independent polynomials of degree k. Let g{t) be the generating 
function for these numbers, 

g{t)=Y,dkt^ (4.5) 

k 

Then this generating function can be calculated as 

In the case at hand the subgroup of the permutation group is actually the permutation group 
of n elements itself. 

The permutation group has a non-trivial action on the space of solutions to the momen- 
tum conservation constraints. This space can be spanned for instance by the set in equation 
(2.6). To illustrate this in the four particle case, one can pick s = {ki + k2)^ and t = {k2 + k^)'^ 
as the independent variables. A permutation of particles one and two then induces 

s^ s t^-s-t (4.7) 

which can be written in matrix notation as 

(4.8) 

It is straightforward to derive the action of all 24 permutations in this case in the same way. 
This yields the Molien series for our problem for four particles as 
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Note that the terms in the Mohen series in equation (4.6) are constant on the different 
conjugacy classes. Hence more efficiently one can sum over only one representative of each 
conjugacy class and multiply by the number of elements in this class 

The Molien series can be used to calculate the generating function for the number of com- 
pletely symmetric, Lorentz invariant polynomials of the external momenta up to momentum 
conservation, ordered by the degree which in this case corresponds to twice the mass dimen- 
sion. More generally, polynomia invariant under any subgroup of the permutation group can 
be obtained this way. Concretely, let the vector space Y be the set of independent Man- 
delstams (2.6) at a fixed multiplicity. The matrix representation of the permutation group 
can be computed for this set. Then equation (4.10) can be used to calculate the generating 
function an the sought-for numbers can then be read from its expansion. 

The explicit expressions for the generating functions get quite complicated quickly. The 
five particle result for instance is 

, , 1 / 24 30 20 



120 V 1 - i^ \-t^ -t^^t 1-t^ -t^-t^ + t^ + t 

20 25 

+ -. ^ r. 7^ 7^ : + 



l-t^ + t^-t^ + t'^-t l-t^ + t^ + 2t^ -2t'^-t 

+ R -. —^ n I (4.11) 

Of course, the above algorithm can easily be implemented in a symbolic computer program 
such as Mathematica. The results for the first 15 degrees of the polynomials for up to and 
including 10 particles are summarized in table 1. Note that read vertically at column i the 
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205 
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4321 
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84927 
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2 
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11 


33 


81 


246 


722 


2222 


6875 


21497 


66299 


202179 



Table 1. Numbers of completely symmetric polynomials up to momentum conservation as a function 
of degree up to and including 10 particles. 

series seems to stabilize at 2i — 2. This we have checked up to 16 particles (disregarding 
Gramm determinant constraints). If this observation turns out to be true, this gives the 
numbers of polynomials up to degree 8 as in table 2. 
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n > 16 



012345 6 7 8 
"l i 2 5 n 34 87 279" 



Table 2. Conjectured asymptotic numbers of completely symmetric polynomials up to momentum 
conservation. 

Basis 

In the four point case a basis for the completely symmetric polynomials is known as displayed 
in equation (4.2). The Molien series can be used to find a candidate for a similar basis 
at higher points. The argument proceeds degree by degree: first, at degree 2 there is only 
one independent symmetric polynomial. Pick one representative of this. If a set of basis 
polynomials is known at degree i, then one can construct all polynomials at degree i + 1 hy 
taking products within the set such that the resulting degree is i + 1. The difference between 
the number of polynomials thus generated and number obtained from the Molien series at 
this degree, if positive, is a lower bound on the number of polynomials which have to be 
added to the basis at degree i + 1. It is a lower bound, as there can be relations between 
the products of basis-polynomials induced by momentum conservation. This algorithm can 
be used to estimate the number of different basis polynomials which are needed to generate 
all other polynomials. The results are listed in table 3 up to seven particles up to degree 18 
polynomials. For four and five particles we have checked up to degree 100 polynomials that 
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247 
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497 
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Table 3. Numbers of needed basis polynomials as a function of degree based on numerology. 

the indicated set generates enough polynomials at higher degree to explain the Molien series 
results. For four particles this just reconfirms the result expressed in equation (4.2). For 
five particles based on the presented numerology there is a candidate basis consisting of nine 
elements. These will be chosen to be 

Ti = (s*i2 + permutations) i G {2, 3 . . . , 9} (4.12) 

as well as 

n = sl^sl^ (4.13) 

We have checked explicitly to degree 9 these nine elements indeed form a basis. Note that this 
does not exclude relations between polynomials generated by this set of nine basis elements 
at higher degree: the first instance of this is at degree 12 where there is one relation between 
the 24 possible polynomials. 
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In conclusion, in this subsection it was argued that for five particles all symmetric poly- 
nomials of the Mandelstams can be expressed in the r basis just presented, e.g. 

9 

h= T. c,„i„...,,ll{nr (4.14) 

il...ig k=l 

If the polynomials / is homogeneous of degree k then 

9 

6h + ^j{ij) = k (4.15) 

should hold^ for the integer valued coefficients ij. This expansion is not necessarily unique 
from degree 13 onwards. 

4.2 Obtaining explicit results for MRV amplitudes at five points 

The stage is now set to compute the MRV amplitudes at five points from the KLT relations 
and then express the result in the basis just constructed. This calculation was alluded to in 
[25]. A generic consequence of the superspace form of the amplitude in equation (4.1) is that 
the pre-factor A multiplying the fermionic delta function follows basically from computing 
one component amplitude. In this case it is easiest to pick the amplitude with four gravitons 
and one holomorphic scalar out of the IIB spectrum^. The scalar is generically a combination 
of axion and dilaton, but the axion part is necessarily zero here (it would arise through KLT 
from open string amplitudes with a single fermion). The needed fermionic integral can be 
computed easily in a four dimensional approach that can be used here since the problem only 
involves five-particle kinematics. The result of this integral is, in spinor helicity notation, 

DC (12)^ [34]^ after assigning helicity -|--|- to legs one and two and to legs three and four. 

This gives 

A{12)'[Uf = A{++, ++, — , — , 0) (4.16) 

The right hand side of this equation may now be calculated by the KLT relations. The scalar 
(j) decomposes in a sum over opposite gluon polarization for the fifth leg. Schematically, this 
reads 

^cl(++,++, , ,<^)~ 

^o(+, +, -, -, -)^o(+, +, -, -, +) + ^o(+, +, -, -, +)^o(+, +, -, -, -) (4.17) 

where the open string amplitudes on the right hand side are color ordered. From the explicit 
expressions for the five point MHV amplitude and its conjugate it is easy to factor out 
the universal (12)^ [34]^ factor. The resulting expression for A is after some rewriting into 
Mandelstam invariants: 

A = ^ fCn + '-^^G,, - ^1^^34 -.23^14 + . 13.24 ^\ ^^^^^^ 

S12S23S34S45S51 V S13S24 S13S24 / 



Note this type of equation is neatly solved in Mathematica by the command "FrobeniusSolve" 
''Calculation at low orders for these and similar amplitudes also appear in [10]. 
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where Gij refers to an entry from the matrix 

G = F^^F = C'^SqC (4.19) 

From this expression it is not obvious this is going to yield a completely symmetric polynomial 
in the external momenta. This is true though in all known instances and is guaranteed by the 
underlying Bose symmetry of the closed string. Given the previous results it is now a simple 
matter to obtain the MRV amplitude up to order 21 by the results obtained above. 

To write the algebraically complex result in a nice format, one can employ the same 
numerical trick as mentioned before for verifying (2.21). This produces the coefficients of 
the field theory expansion of the MRV amplitudes into the 9 element basis. Again, this is 
easily automated. Perhaps more interestingly, the conjecture in equation (2.21) can be used 
to produce amplitudes in the MRV sector to, in principle, arbitrarily high order. In particular 
coefficients of powers of single zeta's are readily isolated. As an example, consider the simplest 
series, 

Ac, = 1 (4.20) 

Aci = -\rs (4.21) 



1 3 1 2 / 5 , 1 



A^3 = — ^To H To - (-^T2T4) + -Tfi - Tl (4.22) 



Higher order results can be obtained but do not offer much insight and are therefore sup- 
pressed. It would be interesting to see if there is more structure to be found exploring these 
ideas. For instance the basis above was chosen for constructional convenience only: it would 
be interesting to see if string theory would prefer one particular basis. 

5 Discussion and conclusion 

In this paper the field theory expansion of superstring theory amplitudes with massless legs 
in a fiat background has been discussed. The main idea which may generalize to higher 
points is to introduce an additional parameter into the integral (the 'z' of the 3F2 function) 
and consider the field theory expansion for generic values of this parameters first. This idea 
certainly generalizes beyond the five point case and it will be interesting to see what further 
consequences may be derived from this. Using the Mellin-Barnes argument given at the start 
of section 3 it is easy to see that above five points more complicated functions than hypergeo- 
metric ones will arise. These are generically boundary values of hypergeometric-like functions 
as discussed in [18]. Concretely it would be interesting to explore a similar differential equa- 
tions based approach to expanding these functions, if only to see what additional structure 
may arise. 

Resolution of the niggling discrepancy at weight 21 between conjecture and calculated 
(j) map is certainly called for. Furthermore, structurally it will be interesting to explore 
further applications of the recursive structure for the field theory expansion of the five point 
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amplitudes. Here it would be most welcome to relate the appearing 6 operator for instance 
to a natural world-sheet quantity. More generally, the structure of the field theory expansion 
encodes deep properties of the string theory and it is certainly worthwhile to explore these 
further, for instance starting with [8]. One direction here is reconstructing the effective action 
from the obtained field theory expansion. 

For MRV amplitudes an extension of arguments used in [10] in the four dimensional case 
to IIB in 10 dimensions would be most welcome to push beyond five legs. This is interesting 
for various reasons, not least because of the known connections between MZVs and modular 
functions [28] which is certain to contain wonderful results waiting to be explored. It is 
expected that the multi-zeta value structure follows most natural in the context of MRV 
amplitudes in the closed IIB string. 

Molien's theorem certainly deserves to be wider known than it is currently in the physics 
community, especially in places where the permutation group makes an appearance. From the 
discussion above it should be obvious how to count and analyze much more general problems. 
Finally, it is hoped that the explicit results obtained in this article and included in the archive 
submission can be useful for others for applications beyond the ones mentioned here. 
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